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We study the relation between energy and entanglement in an entanglement transfer problem. We 
first analyze the general setup of two entangled qubits ( "a" and "6" ) exchanging this entanglement 
with two other independent qubits ["A" and "B"). Qubit "a" ("6") interacts with qubit 11 A" ("B") 
\Q . via a spin exchange-like unitary evolution. A physical realization of this scenario could be the 

problem of two-level atoms transferring entanglement to resonant cavities via independent Jaynes- 
Cummings interactions. We study the dynamics of entanglement and energy for the second pair of 
qubits (tracing out the originally entangled ones) and show that these quantities are closely related. 
For example, the allowed quantum states occupy a restricted area in a phase diagram entanglement 
vs. energy. Moreover the curve which bounds this area is exactly the one followed if both interactions 
are equal and the entire four qubit system is isolated. We also consider the case when the target 
pair of qubits is subjected to losses and can spontaneously decay. 

On 

£Nj . PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Yz 



> 
On 

Cn ■ I- INTRODUCTION 

oo 

Entanglement is one of the most studied topics at present. The large interest for this issue relies mainly on the fact 
that entangled systems can be used to perform some tasks more efficiently than classical objects It is then natural 
■ to look for a good understanding of this resource not only from a purely mathematical point of view, i.e., formalizing 
""^5 | the theory of entanglement, but also from a more practical approach, i.e., studying its role and manifestations in 
realistic systems. For example, recent works have been able to connect entanglement to thermodynamical properties 
of macroscopic physical systems |2j. In a distinct venue, other works study physical manifestations of quantum 
correlations by suitably choosing particular purity and entanglement quantifiers and restricting allowed quantum 
states according to these quantities Q- In these studies, concepts like maximally entangled mixed states (MEMS) 
are discussed. A very recent study also adds energy to entanglement and purity as a third parameter to characterize 
certain quantum states 0. In particular, the authors discuss the physically allowed states according to the possible 
values of entanglement, purity and energy for a system composed of two qubits or two gaussian states, and also study 
the entanglement transfer between them. 

In the present manuscript, we study the connection between entanglement and energy that appears naturally in 
a swapping process involving two systems of two qubits. In the model investigated, we consider a simple form of 
interaction between two pairs of qubits labeled as a A and bB. The system ab is prepared in an entangled state while 
the pair AB is prepared in a factorable state. We analyze the dynamical relations between energy and entanglement 
of qubits AB when exchanging energy and coherence with qubits ab. In particular, for any given time t, we calculate 
the full quantum state of qubits abAB and then we trace out qubits ab to calculate energy and entanglement of the 
remaining pair AB. We show that this dynamics yields paths in an entanglement-energy diagram, and that these 
paths are contained in a very restricted region. Moreover, we identify the frontiers of this region from the general form 
of the density operator that represents the state of the subsystem AB. We also propose a physical system to realize 
such entanglement transfer and investigate how the dynamics of the entanglement swapping is modified if the AB 
system is open and allowed to dissipate energy to an external reservoir. In some sense, this work is complementary to 
the sequence Q in which the authors study the problem of entanglement transfer from continuous- variable entangled 
states to qubits, although in those works the authors do not pay particular attention to the relation between Energy 
and Entanglement. 

The paper is organized as follows. In section [n] we introduce the general physical system that we will study and 
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the basic setup from which we will approach it. We also define the quantities that will be analyzed throughout the 
paper and finally we discuss the dynamics of this system. Section Ifffl is devoted to study the entanglement and the 
energy of system AB under a particular unitary evolution. We then propose a physical implementation for the studied 
Hamiltonian, and generalize the time evolution in section llVl by considering the problem of a dissipative, non-unitary 
evolution. In section[3we conclude by reviewing the main points we have discussed and suggesting possible extensions 
of this study. 



II. PHYSICAL SCENARIO 



we are interested in. Suppose a system of four qubits a, b, A, and B interacting 



Let us start by describing the system 
via a spin-exchange like Hamiltonian: 
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For each qubit, the relevant Pauli operators are defined by 

a z = |1> <lj — ]0> <0| , (3a) 
a+ = |1}(0|, (3b) 
a. = |0)(1|, (3c) 

and the interaction operators like ata^, for example, can be viewed as annihilating an excitation of subsystem a and 
creating an excitation in subsystem A. The constants g a A and gbB give the strength of the interaction between these 
subsystems. One important feature in understanding such Hamiltonians is that the total number of excitations is a 
conserved quantity. The eigenvectors of (|2a|l (similarly to (|2b|^ are given by: |00) aj4 , with eigenvalue E^q = —Hut, 
\ll) aA , with eigenvalue Eff = fuv, and \A±) aA = (|01) ± \10))/s/2, with eigenvalue E± A = ±hg aA /2, where uj = 
(uja +uja)/2- 

As the initial state, let us suppose that the entire abAB system is prepared in the form: 

|<^ = O)) = |V(0)) o6 ®|OO> AB . (4) 

where \ip(0)) ab — sin 9 |01) +cos# |10), which means that subsystem AB is prepared in its ground state and subsystem 
ab is usually prepared in some entangled state with one excitation (except if 9 = n^, n G Z, when the state is 
factorable). Note that this initial state is pure and it is chosen so that the bipartition ab <g) AB does not present 
any initial entanglement. From now on, we will study the time evolution of this initial state when subjected to 
Hamiltonian JQl for different coupling constants g a A and gbB- We will concentrate our analysis in the subsystem AB 
by tracing out the degrees of freedom of systems a and b. Another simplifying assumption we made is to consider the 
complete resonance condition uj a — uja = ujb = ujb = lj- 

A special case of this dynamics happens when g a A = gbB = g, in which case state (0J evolves into state 

W)) = cos(gt) \^) ab (g) |00) AB - i sm(gt) \00) ab ® \^} AB . (5) 

Note that in this simple case, for t = n^, with n odd, the subsystems exchange their states, the entanglement initially 
present in subsystem ab is completely transferred to subsystem AB and with respect to the bipartition ab Cg> AB the 
state becomes again separable. However, for t ^ n^, the whole system is entangled (as long as 9 ^ n^) and subsystem 
AB will be in some mixed state. 

In a more general situation (different coupling constants), state (@J will evolve into: 



|0(i)) = cos9[cos{g aA t)\1000}-ism(g aA t)\00l0)} 
+ sm9[cos{g bB t) |0100) - isin( 56S f) |0001)]. 



(G) 
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Note that for generic times t, state © presents, again, multipartite entanglement among all its individual components 
(a, b, A and B). Studying this multipartite entanglement may also prove intriguing and enlightening. However this 
is not the purpose of this manuscript where, as mentioned above, we will concentrate our analysis in the subsystem 
AB. 

Our goal is to investigate the relation between energy and entanglement in subsystem AB as a function of coupling 
constants and time. In order to study entanglement we will use the negativity [N) which can be defined for two qubits 
as two times the modulus of the negative eigenvalue of the partial transposition of the state p, p TA [|J, if it exists. 
For short: 

N(p) = 2max{0,-A mm }, (7) 

where X m i n is the lowest eigenvalue of p TA . Our choice is motivated by the facts that the Negativity is easy to 
calculate and provides full entanglement information for a two-qubit system. For the energy of subsystem AB, U, we 
will consider the mean value of the relevant restriction of the free Hamiltonian: 

U = Tr{pH AB } 1 (8) 

where 

Has = ^ {at + of) . (9) 



III. ENTANGLEMENT AND ENERGY 



After tracing out the degrees of freedom of systems a and b in the global quantum state © , the reduced state for 
the pair AB is described by: 
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(10) 



where a + b + c= 1 (from the normalization of pab), with a, b, c and d given by the following functions of the coupling 
constants and time: 



a = cos 2 9cos 2 (g a At) + sin 2 9 cos 2 (gi, B t) , 

b = sin 2 9 sin 2 (gt,Bt) , 

c = cos 2 9sm 2 (g aA t), 

d = cos 9 sin 9 s\ii{g a At) sin^&sf). 



Following Eq. JHJ, the energy of state (fTU|) is: 
which, by means of Eqs. (JUJ, becomes: 



U = 



U 



cos 2 9 cos 2 (g a At) — sin 2 9 cos 2 (gbBt)- 



(11a) 
(lib) 
(11c) 
(lid) 



(12) 
(13) 



Note that — 1 < U < 0, which means that there is at most one excitation on the AB system. This is expected since 
the chosen initial state contains only one excitation for the entire abAB system and this set of qubits is isolated, i.e. 
it cannot be excited by external sources. A simple calculation gives for the entanglement (negativity) of state IjlUfl 



N = vV + 4rf 2 



(14) 



A. Entanglement and energy versus time 

In Fig. we have plotted the temporal behavior of N and U for several values of the coupling constants g a A and 
g^B- The entanglement transferring process can be followed in those pictures. The simplest one is for equal coupling 
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N,m=n=1 
N,m=1,n=2 
U,m=n=1 
U,m=1,n=2 




N,m=1,n=3 
N,m=3,n=5 
U,m=1,n=3 
U,m=3,n=5 



1i 




N,m=1,n=sqrt(2) 

N,m=1 ,n=sqrt (3) 

U,m=1 ,n=sqrt 2) 

U,m=1,n=sqrt(3) 

FIG. 1: (Color online) Energy U and Negativity TV of the state pab given by Eq. IIUH versus time. The initial state is given by 
Eq. Q with 9 — 7r/4, which means a maximally entangled pair ab is initially present, and its entanglement can be transferred 
to the pair AB. Several values of m and n were used in the relation mg a A = ngbB- 

constants g a A — gbB, in which the state \ip) is cyclically "bouncing" between the two pairs of qubits. The chosen ratios 
between coupling constants indicate a very important behavior of the system: the complete entanglement transferring 
process can only happen iimg a A — ng^B for m and n odd integers. Eq. Ijllal) supports this conclusion, since both g a At 
and gbst must be odd multiples of ? simultaneously (remember that cos 2 6 and sin 2 9 are positive numbers in order to 
the state \ip (9)} to be entangled). The physical picture is that each pair aA and bB oscillates inside duplets 1 10} and 
1 01) . The situation is analogous to two classical harmonic oscillators with distinct frequencies starting from a common 
extremal point. The implied relation is necessary for them to meet within an odd number of half oscillations, which 
is the condition for a complete transfer of state. To insist in this point, note that for g a A = ZgtB, at time t = 
the pair aA (or, more precisely, its analogous oscillator) has suffered a full oscillation, but the pair bB (resp. its 
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analogous oscillator) has undergone half an oscillation, so one could found entanglement between a and B, but no 
entanglement can be found between any other pair of qubits, including the studied pair AB. 

For other values of 9 we obtain similar pictures, with the only important difference in the entanglement scale, since 
no maximally entangled pair will be formed. 

B. Entanglement versus energy 

In this subsection we use time as a parameter to draw graphics on an entanglement versus energy diagram. As we 
will show, the paths followed in this phase diagram exhibit interesting patterns. 




FIG. 2: (Color online) Negativity (N) versus Energy (U) for pAB given by Eq. 1 1 (H with 9 — n/4, i.e. the initial ab state is 
maximally entangled. Parameter relation g a A = ng^B with n = 1 (left upper panel), n = 2 (right upper panel), n — 7 (left 
lower panel), n — 53 (right lower panel). 

The graphics for 9 = ir/4 (Fig.|2J| and 9 = tt/3 (Fig.[5J are qualitatively different. However, it can be seen through 
Fig. 0] that, independent of the available initial entanglement (given by the value of 9), the accessible region in the 
parameter space N x U is bounded by an upper curve. This bound can be explained in the following way: using the 
fact that U = — a in Eq. (|14|l we have 

N 2 - 2NU = 4\d\ 2 . (15) 

However, when g a A — gbB and 9 = 7r/4 we have b = c = d, in which case iV 2 — 2NU = Ab 2 . Using the normalization 
condition a + b + c = 1, we get 

(l + U) 2 =M 2 . (16) 

Therefore, in the ideal situation of equal coupling g a A = gbB and maximally entangled initial ab state (9 = 7r/4), we 
can write 

N 2 — 2NU = (1 + U) 2 . (17) 
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FIG. 3: (Color online) Negativity (N) versus Energy (U) for pAB given by Eq. I| !()> with 9 = 7r/3, i.e. the initial ab is partially 
entangled. Parameter relation g a A = ngtB with n — 1 (left upper panel), n — 2 (right upper panel), n — 7 (Left lower panel), 
n = 53 (right lower panel). 



As we will see now this equation is exactly the one that limits the phase-space for quantum states in this problem. 
The normalization condition yields b + c = 1 + U that allows us to obtain 

4bc={U + l) 2 -(b~c) 2 , 

which implies 

4bc<(U + l) 2 . (18) 

At the same time, the condition for matrix IjlUI) to be considered a true density matrix is that its eigenvalues are all 
positive, which is reached if and only if \d\ 2 < be. Therefore, we can conclude that 

4|d| 2 < 46c < (U + l) 2 . (19) 

and, from Eq. JTBJ, we find 

N 2 -2NU< (U + l) 2 . (20) 

As we saw in Eq. (II 711 the equality is reached for g a A = gbB and 9 = n/4. So Eq. I|17f) bounds the region that density 
matrices of the form (|10|) can occupy in the diagram N x U . 



IV. OPEN SYSTEM 



Up to this moment we have considered any two pairs of qubits. Now we will adhere to one specific physical 
realization, namely: two atoms resonantly coupled to two independent cavity modes. If the cavities are initially in the 
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ground state (vacuum, no photon) and the usual approximations are valid Q, the Jaynes-Cummings Hamiltonian 

Hjc = ^-v z + (a^a + ij + Ml, (aV- + acr + ) (21) 

essentially reduces to the form , with the lower case qubit representing the two- level atom and the capital one, the 
first two energy levels of the field mode. Hence, the situation here studied models an experiment where previously 
entangled atoms transfer such entanglement to independent cavity modes. One nice point when considering this 
particular situation of atoms transferring entanglement and energy to resonant cavities is that those systems can have 
very different dissipation times. In fact, atomic levels can be selected so that their dissipation time scale is much 
larger than those of typical resonant cavities. In this case, one can ask what happens if the system that receives 
the energy and the entanglement dissipates it to an external reservoir. In order to answer this question, the unitary 
analysis considered up to now has to be abandoned and we must change from a Hamiltonian approach to a master 
equation one. 

We will consider that the qubits AB, now represented by the cavity field modes, are in contact with independent 
reservoirs and interact with them. Since atomic lifetimes (for the atomic transitions used in cavity QED experiments) 
are usually much greater than cavity decay times, we will not couple the lower case qubits to any external device. 

To address this problem we consider the time evolution of the global system described by the master equation in 
the Lindblad form 

J t PaA bB {t) =^[H, PaAbB (t)] + l-Y,([^PaAbB(t),V^] + [% PaAbB (i) V?] ) , (22) 

i 

where H is given by and the operators V and V} describe the effects of the coupling to the reservoirs. For 
simplicity, we will model only the dissipation of energy in the cavities coupled to null temperature reservoirs, which 
can be done using 

Vi = y / 2tiK A a^ <g) I abB (23a) 
% = ^/2hn B I aAb ® ct b (23b) 
Vi = , V i > 2 (23c) 

The constants n A and k b are directly given by the decay rates of each cavity mode. 

Starting from the special case of the initial state Q, given by = j (a Bell state for the donor pair of qubits), after 
a somewhat lengthy calculation we find that the state of both cavities can still be written in the form of Eq. JTDJl, but 
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now with the matrix elements given by 




■ sin 



\ 2 9aA „:_^aA t y- KAt /2 
gbB . f^bB 



-1 2 



V2 
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sin 



bB 



V% g a A . (QaA 
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llbB V 2 



with the definitions 



bB 



45, 
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(24a) 

(24b) 
(24c) 
(24d) 

(25a) 
(25b) 



Since the state of the system AB is still described by density matrices of the form we expect the existence of the 
same bounds in the energy-time diagram H17|l . Also note that energy and negativity are still respectively described by 
equations of the forms l|12f) and {H}. In Fig. [5]we have plotted some curves of negativity and energy vs. time, for the 
system AB considering the temporal evolution of the matrix elements given by Eq. I|24[) . Note that the graphics are 
qualitatively similar to the ones displayed in Fig.nf° r unitary evolutions. However, as expected, both entanglement 
and energy decay exponentially to their lowest values as a function of time. This can be understood from the fact 
that the environment drives the system exponentially to the state |00) (see the matrix elements in Eq. (|24H - the 
element (|24a|) goes to unity while all the others go to zero) , which has no entanglement and has the minimum energy 
value U = — 1. 

We have also plotted the negativity vs. energy for the non- unitary case for different values of g a A, 9bB, i^aA and 
Kj,b- This is displayed in Fig. [fj] As commented before the followed paths are still bounded by the same limits of 
the unitary case. However, as the dissipative mechanisms get stronger (i.e. coefficients ka and kb get closer to the 
coupling constants g a A and gbB), less and less entanglement is transferred to subsystem pab- 

A different phenomenology can be anticipated for the case of non-null temperature. Since thermal photons can 
now be captured by both cavities, the state |11) will be populated and the form IjlUfl will not be valid anymore. One 
consequence is that one can expect the phenomenon known as entanglement sudden death since it will not be 
necessary to nullify the elements d in order to have a positive partial transpose, and hence no entanglement. 



V. DISCUSSIONS 



In this paper we have addressed the problem of entanglement and energy transfer between pairs of qubits. We 
considered the particular example of two atoms interacting with two cavities in the Jaynes-Cummings model. This 
evolution can be seen as a state transferring process and, for specific coupling constants, a dynamical entanglement 
swapping. If the atoms are initially in an entangled state, this entanglement is fully or partially transferred to the 
cavities depending on coupling constants and time. This entanglement swapping process is accompanied by an energy 
transfer as well, and we have shown that entanglement and energy in the cavities system are strictly related. 

To clarify this relation we studied these quantities in various scenarios. First we considered the whole system as 
isolated, and investigated its time evolution for several coupling constants and different initial atomic entanglement. 
In each case, we traced out the atoms (as we now refer to the lower case qubits) , we drew an entanglement vs. energy 
phase-diagram for the cavity modes and we found an upper-limit for all the possible paths in these diagrams. This 
bound corresponds to maximally entangled atoms transferring its entanglement and energy to independent cavities 
at exactly the same rates. 

We also considered the possibility of dissipation in the cavities. In this case, while the atoms are transferring excita- 
tions and entanglement to the cavities, some energy is lost to the environment. The cavities state goes asymptotically 
to state 1 00). In the dissipative regime, the evolution of entanglement and energy of the cavities state exhibits the 
same characteristics pointed out for the unitary case, i.e., the paths followed in the entanglement-energy diagram are 
limited to a restricted region whose frontier is identified by the trajectory described when the couplings are identical 
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FIG. 5: (Color online) Energy U (red) and Negativity N (blue) of the cavity modes state (1I0> with matrix elements given 
by Eq. {3]J versus time. The initial state of the system is given by @ with 9 — n/4, i.e., a maximally entangled ab state. 
Parameter relations: n a A ~ KbB ~ 0.1g a A and gbB/g a A = n. Left-above: n = 1. Right-above: n = 2. Left-below: n = 3. 
Right-below: n — \[2. 



and the initial entanglement of the atomic state is maximum. However, as expected, neither entanglement nor en- 
ergy can be fully transferred unless the dissipative times are much larger than the inverse Rabi frequencies involved. 
Fortunately, this regime is usually achieved in cavity QED experiments. 

We only analyzed the entanglement between qubits AB (the modes, in the physical realization proposed). However, 
most of the time the whole system presents multipartite entanglement which may provide interesting new results if 
further studied. Other important continuations of this work include the treatment of non-completely resonant systems 
(e.g.uja = oja 7^ tob = which corresponds to two distinct atoms resonantly coupled to cavity modes, as well as 
the case of dispersive coupling) and also other couplings to reservoirs, like including temperature in the scenario here 
presented and also considering spontaneous decay for the atoms. 



Acknowledgments 

The authors recognize fruitful discussions with M.C. Nemes. Financial support from CNPq and PRPq-UFMG is 
acknowledged. This work is part of the Milleniun Institute for Quantum Information project (CNPq). 



[1] M. A. Nielsen, I. L. Chuang, Quantum Computation and Quantum Information. (Cambridge Univ. Press, Cambridge, 2000). 

[2] X. Wang and P. Zanardi, Phys. Lett. A 301, 1 (2002); C. Brukner and V. Vedral, |quant^ph /0406040 M. R. Dowling, A. 
C. Doherty, and S. D. Bartlett, Phys. Rev. A 70, 062113 (2004); O. Giihne, G. Toth, and H. Briegel, New J. Phys 7, 229 
(2005); O. Giihne and G. Toth, quant-ph/0510186 I. Bose and A. Tribedi, Phys. Rev. A 72, 022314 (2005); G. Toth, Phys. 
Rev. A 71, 01030RR) (2005); M. Wiesniak, V. Vedral, and C. Brukner, New J. Phys. 7, 258 (2005); C. Brukner, V. Vedral, 
and A. Zeilinger, Phys. Rev. A 73, 012110 (2006); T. Vertesi and E. Bene, Phys. Rev. B 73, 134404 (2006); 



10 




FIG. 6: (Color online) Negativity (N) vs. Energy (U) of the cavity modes state 11 1 01 with matrix elements given by Eq. JSJ for 
fixed decay rates and different couplings. The initial state of the system ab is given by (|1J with 6 — 7r/4. Parameter relations: 
K a A = KbB ~ 0.1g a A and g a A = ngbB with n = 1 (left upper panel), n — 2 (right upper panel), n = 7 (Left lower panel), n — 53 
(right lower panel). 



[3] S. Ishizaka and T. Hiroshima, Phys. Rev. A 62, 022310 (2000); F. Verstraete, K. Audenaert, T. De Bie, and B. de Moor, 
Phys. Rev. A 64, 012316 (2001); T. -C. Wei et al, Phys. Rev. A 67, 022110 (2003); M. Ziman and V. Buzek, Phys. Rev. 
A 72, 052325 (2005). 

[4] D. McHugh, M. Ziman, and V. Buzek, e-print quant-ph/0607012 

[5] W. Son et al, J. Mod. Opt. 49, 1739 (2002); M. Paternostro et al, Phys. Rev. A 70, 022320 (2004); A. Serafini et al, Phys. 
Rev. A 73, 022312 (2006). 

[6] J. Eisert, PhD thesis University of Potsdam, January 2001; G. Vidal and R.F. Werner, Phys. Rev. A 65, 032314 (2002); K. 

Audenaert, M. B. Plenio and J. Eisert, Phys. Rev. Lett. 90, 027901 (2003). 
[7] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, England, 1995). 
[8] G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

[9] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 (2004); L. Jakobczyk and A. Jamroz, Phys. Lett A 333, 35 (2004); 
M. Yonac, T. Yu, J. H. Eberly, J. Phys. B: At. Mol. Opt. Phys. 39 S621-S625 (2006). 



